DNA methylation of imprint control regions associated with Alzheimer’s disease in non-Hispanic Blacks and non-Hispanic Whites

Alzheimer’s disease (AD) prevalence is twice as high in non-Hispanic Blacks (NHBs) as in non-Hispanic Whites (NHWs). The objective of this study was to determine whether aberrant methylation at imprint control regions (ICRs) is associated with AD. Differentially methylated regions (DMRs) were bioinformatically identified from whole-genome bisulfite sequenced DNA derived from brain tissue of 9 AD (5 NHBs and 4 NHWs) and 8 controls (4 NHBs and 4 NHWs). We identified DMRs located within 120 regions defined as candidate ICRs in the human imprintome (https://genome.ucsc.edu/s/imprintome/hg38.AD.Brain_track). Eighty-one ICRs were differentially methylated in NHB-AD, and 27 ICRs were differentially methylated in NHW-AD, with two regions common to both populations that are proximal to the inflammasome gene, NLRP1, and a known imprinted gene, MEST/MESTIT1. These findings indicate that early developmental alterations in DNA methylation of regions regulating genomic imprinting may contribute to AD risk and that this epigenetic risk differs between NHBs and NHWs. Supplementary Information The online version contains supplementary material available at 10.1186/s13148-024-01672-4.


Introduction
More than six million Americans are affected by Alzheimer's disease (AD) [1], which is the most common form of dementia (60% to 80% of cases) [2,3], and it is now the sixth leading cause of death in the USA [3].Additionally, AD places a tremendous burden not only on the patients, but also on the caregivers and the healthcare system.AD is a disease of cognitive changes, but also increases susceptibility to multiple comorbidities, including pneumonia, femur fractures, and increased mortality risk [3].This neurodegenerative disease is associated with cell death and atrophy involving various brain regions, which progress along anatomically connected networks, starting in the entorhinal cortex and medial temporal lobes, and extending into the neocortex over time.Although it is now accepted that this neuropathology is characterized by the aggregation of extracellular amyloid-beta (Aβ) plaques and intracellular neurofibrillary tangles (NFTs) composed of the hyperphosphorylated tau protein [4], the etiologic factors contributing to AD are still largely unknown.Established risk factors for AD, including advanced age, familial history, genetics, history of head trauma, and cardiovascular diseases do not fully explain the formation of Aβ plaques and NFTs [1].
However, up to 95% of the disease is estimated to be sporadic [16,17].Such cases share common neuropathological endpoints with familial AD, including Aβ plaque accumulation, NFTs, synaptic loss, excess inflammation, oxidative damage, and neuronal death [18].While genetic factors, such as APOE variants, appear to influence sporadic AD through intricate interplay with each other and environmental influences, it is important to note that they are neither necessary nor sufficient for the development of AD [17].This has led to the emerging hypothesis that environmental stressors accumulated over the life course contribute to the later development and progression of AD [19].
Globally, there is a slight geographic variation in AD prevalence [20].Among Americans aged 65 years and older, the risk of AD is twofold higher in non-Hispanic Blacks (NHBs) and 1.5-fold higher in Hispanics compared to non-Hispanic Whites (NHWs) [1].The causes for this disparate outcome are presently unclear.Although a single nucleotide polymorphism (rs115550680) of ABCA7, which regulates lipid transport [21], was recently associated with late-onset AD in a NHB population but not NHW or Hispanics [22], this known genetic variant does not fully explain the higher disease burden in NHBs.A plausible hypothesis to account for the elevated rates of disease in NHB and Hispanic populations is that environmental or life course stressors, such as migration and segregation [23], inadequate medical surveillance, and living in polluted environments, are more common in these populations and result in functional and enduring alterations in the epigenome [24,25].Thus, profiling epigenetic marks that link established risk factors to AD holds promise for early detection, and for identifying novel mechanistic pathways contributing to AD.
Epigenetic dysregulation, which can cause alterations in gene expression in response to environmental stressors, may cause long-term changes in molecular pathways contributing to AD.Indeed, it was recently documented that the average 5-methylcytosine level is decreased in the entorhinal cortex of individuals with AD compared to that in controls [26].DNMT1, a critical factor in the maintenance of DNA methylation, and MeCP1/ MBD2, components of the methylation complex, are also significantly decreased in the entorhinal cortex of AD individuals compared to controls [27].Furthermore, brain-derived neurotrophic factor (BDNF), which functions in cortical neuron maintenance, has increased promoter CpG methylation in both AD brain tissue and blood [28], in support of similarities between methylation patterns in the blood and brain tissues of AD and other dementia patients [29,30].Nevertheless, the interpretation of these data is complicated by lack of replication, and the possibility that methylation levels may change between tissue/cell type and throughout life.
In epidemiological studies, DNA methylation is frequently identified in accessible peripheral blood [31][32][33][34], but because these epigenetic marks can differ between tissues and cell types, they do not always correlate with those from inaccessible cells of affected brain regions.Moreover, because epigenetic marks respond to various environmental cues throughout life, causality is difficult to discern.One exception to these issues is the repertoire of methylation marks controlling genomic imprinting-the human imprintome-which epigenetically regulates the expression of imprinted genes crucial to tissue development during the intrauterine period [35].After fertilization, DNA methylation of imprint control regions (ICRs) in the primordial germ cells (PGCs) undergoes complete erasure, and after sex determination, these regions are remethylated in a sex-dependent manner.This time frame is a window of high susceptibility to epigenetic perturbations due to environmental exposures and stressors that can alter the methylation of these ICRs in PGCs [36,37].PGCs with aberrant methylation can then transfer altered gene expression to the next generation and because of mitotic heritability, this aberrant methylation is conserved in all cell types and tissues in the offspring resulting in altered health effects over the life course [38], and increased susceptibility in adulthood to diseases, such as AD.
Thus, the complete mapping of the human imprintome that is susceptible to environmentally influenced alterations is key to understanding the non-genetic factors in complex diseases [35].It is also important to distinguish between the non-imprinted epigenetic-controlled regulatory sites, and the ICRs involved in regulating the parental-dependent expression of imprinted genes.DNA methylation patterns at non-imprinted sites are cell typespecific and can be responsive to environmental cues throughout life.In contrast, the inherited ICRs, or the somatic ICRs that occur at the stem cell stage of embryonic development, should have the same stable methylation status across all tissues throughout life, including peripheral blood cells and the brain [39].As changes in the brain are likely to start decades before clinical symptoms of AD appear [40,41], the purpose of this study was to use genome-wide approaches to comprehensively identify dysregulated ICRs [35] associated with AD that trace their origin to adverse events in early development.The consistency of imprinted methylation marks across tissues and cell types makes them attractive as early epigenetic biomarkers for AD obtainable from accessible tissues.

Patient and sample characteristics
Characteristics of the 17 NHW and NHB individuals who donated the brain samples used in this study are shown in Additional file 1: Table S1.AD cadavers ranged in age from 63 years to > 89 years (median 84 years), and controls ranged in age from 56 to 87 years (median 74 years).All nine AD samples were obtained from the temporal cortex.Five of eight control brain samples were obtained from the temporal cortex while three were obtained from the cerebellum.The diagnosis of AD was made postmortem through a comprehensive neuropathologic evaluation (Additional file 1: Table S1).The brain tissues used in this study were obtained from the Joseph and Kathleen Bryan Brain Bank at Duke University, which has historical significance, as it contained brain samples that were instrumental in the discovery of the association between APOE-ε4 and late onset and sporadic Alzheimer's disease [5,6].

Association of ICRs with Alzheimer's disease
Bisulfite conversion rate of AD and control samples taken from individuals of NHBs and NHWs showed a > 97% bisulfite conversion in all sample groups.Quality controls revealed no sequence duplication bias (Additional file 1: Fig. S1a), and sequence coverage between 15X and 36X (Additional file 1: Fig. S1b, c).We identified the CpG methylation ratio from replicate bam files for AD samples and controls.Later, differentially methylated regions (DMRs) were called using model-based analysis of bisulfite sequencing data criteria (MOABS, version 1.3.8.7) [42].We performed three different analyses to identify AD-associated DMRs (Additional file 2: Table S2) using a 10% differential methylation threshold, a minimum read depth ≥ 7, and a maximum distance between consecutive CpG's ≤ 300, consistent with MOABS [42].The resulting set of AD-related DMRs were analyzed against the 1488 candidate ICRs reported by our group [43], resulting in the identification of 120 candidate ICRs, including four of the 25 confirmed ICRs, that exhibit differential methylation in AD patients compared to controls (Fig. 1a, Table 1).Stratified by group, I. 40 (33.3%)differentially methylated ICRs are observed between all AD samples (n = 9) and controls (n = 8), II.81 (67.5%) are observed between NHB-AD cases (n = 5) and controls (n = 4), and III.27 (22.5%)are observed between NHW-AD cases (n = 4) and controls (n = 4).Interestingly, our results indicate that NHBs exhibit a threefold increase relative to NHWs in AD-associated differential methylation of regions postulated to be ICRs (Fig. 1a, Table 1).Alignment of AD-related DMRs and candidate ICRs can be accessed at https:// genome.ucsc.edu/s/ impri ntome/ hg38.AD.Brain_ track.
When we constrained the analysis of the 1488 ICRs, specifically to the 332 for which gametic methylation patterns were available [35], we did not observe differences that would affect the proportional differences between AD and controls in both NHBs and NHWs.For example, we identified 37 of 332 candidate ICRs associated with AD, with 31 ICRs found in NHBs and two ICRs found in NHWs (Additional file 1: Table S3).Notably, we again observed a common ICR, ICR_481 (MEST/ MESTIT1), which exhibited differential methylation in AD brains compared to controls.Despite observing a reduced number of DMRs associated with AD when constraining the analysis to the 332 ICRs for which there is gametic data, the demonstration of parental origin of methylation strengthens our confidence in the 37 identified differentially methylated regions as robust  [40], NHB [81] and NHW [27] when DNA methylation differed by ≥ 10% between AD cases and controls.c Venn diagram of ICRs from ALL [10], NHB [32] and NHW [10] when DNA methylation differed by ≥ 15% between AD cases and controls.Created with BioRender.comcandidates for regulatory regions in the development and progression of AD.We identified 106 genes in closest proximity to the 120 AD-associated ICRs (i.e., between 0 and 186,698 bp; average 11,270 bp) (Table 1).This is well within the range of known imprinted domains.For example, the H19/ IGF2 and KCNQ1 imprinted domain is about 1.4 Mb long; that of MEST/MESTIT1 is around 4.0 Mb; and that of the NPP5F_v2 is nearly 8.6 Mb DNA [55].
When stratified by race/ethnicity, there were 85 NHB-AD associated ICRs and 26 were linked to NHW-AD associated ICRs.Network analysis conducted separately on the 85 ICRs in NHBs and 26 ICR in NHWs using ingenuity pathway analysis (IPA) unveiled shared functions such as cell signaling, cellular development, embryonic development, and organ development in both NHBs and NHWs (Additional file 1: Table S4, S5).Interestingly, two pathways, namely white adipose tissue browning and gap junction signaling, were also identified as common features in both NHBs and NHWs.On the other hand, the netrin signaling pathway, known to regulate axonal growth, was found only in NHBs.

Discussion
Although it is established that mutations in APOE, APP, PSEN1/2, and BACE1 contribute to AD risk [77,78], and that the APOE-ε4 allele affects cognitive function [77,79], known genetic variation alone explains only a small proportion of AD.Evidence in the last decade supports that epigenetics may contribute substantially to altered gene function and disease development.The highly conserved The candidate ICRs were previously defined by having 5 or more consecutive CpGs with methylation levels of 50% ± 15% (green dots) for tissues in all three germ layers (i.e.brain, kidney, and liver); methylation levels for sperm and oocytes are also shown (i.e.≥ 90% methylation-yellow dots and ≤ 10% methylation-blue dots) [35] and stable methylation pattern of ICRs makes them particularly valuable in the study of diseases like AD, which do not manifest until adulthood or advanced age.After fertilization and remodeling of methylation status during the intrauterine period, ICRs normally maintain the same methylation status in all cells and tissues, including in blood and brain tissue, throughout life.For this reason, early changes in the methylation of ICRs could potentially serve as susceptibility biomarkers for disease risk, and they could be measured at any time during an individual's life.To elucidate the AD associated candidate ICRs in brain tissue that have potential regulatory functions, we determined the methylation pattern of the genome by WGBS and identified the differentially methylated regions in AD brain samples compared to that in controls.We identified 120 candidate ICRs with altered methylation levels in patients with AD.
We next determined whether the patterns of ADrelated methylation in candidate ICRs differ by racial/ ethnic group.A threefold difference in the number of AD-related ICRs was found in the brain samples of NHBs (67.5%) when compared to NHWs (22.5%).This finding is consistent with the postulate that environmentally responsive epigenetic differences in the methylation of ICRs could contribute to the racial/ethnic disparities observed in AD between NHBs and NHWs [1].
Remarkably, one of the two common ICRs identified in NHBs and NHWs was ICR_481 (MEST/MESTIT1).MEST/ PEG1, a paternally expressed imprinted gene, was shown to regulate neuronal migration in development of neocortex [80,81], block neuron differentiation when knocked out, and inhibit Wnt signaling when expressed [50].This is functionally significant, as Wnt signalling has been reported to be associated with age related neurodegenerative diseases including AD [80,82].As MEST expression is reduced by promoter hypermethylation, this would result in activation of Wnt signaling in brain tissues of AD patients, potentially facilitating the progression of AD [80].
The second ICR identified in NHBs and NHWs was ICR_987.It is closest to the gene NLRP1, a component of the inflammasome complex that triggers an immunostimulatory form of cell death called pyroptosis, which is activated in neuronal cells in response to amyloid-β (Aβ) aggregates [53,83].Inflammasomes are multiprotein complexes that are assembled in response to a cellular stressor, including infection, and lead to caspase activation [84].They also are associated with neurodegenerative diseases such as AD [53], consistent with the hypothesis that neuroinflammation contributes substantially to neurodegeneration [53].Studies in murine AD models indicate that the Nlrp1 inflammasome is indeed upregulated, and neuronal death is observed, leading to cognitive decline [53].Kaushal et al. [54] reported a 25to 30-fold higher number of NLRP1-expressing neurons in AD brains compared to control brains.The existence of significantly increased methylation of the candidate ICR_987 (NLRP1) in both NHBs and NHWs supports the neuroinflammation hypothesis of AD formation, one of the most studied mechanisms in AD pathogenesis.Furthermore, there are four single nucleotide polymorphisms (SNPs) (rs2137722, rs11657747, rs34733791 and rs11651595) in the NLRP1 region that have been reported as significantly associated with AD [85].The major alleles for SNPs rs2137722 and rs11657747 are Gs in CpG sites, such that the minor alleles abolish potential methylation sites.Interestingly, the minor A allele for rs2137722, which would block methylation, appears to provide a protective effect against AD [85].This functional importance of NLRP1 in the development of AD further supports the potential regulatory importance of the proximal differentially methylated ICR_987.
We conducted a comparative analysis of data from three previous epigenome studies: Zhang 2023) in blood [88].We aimed to identify overlaps between previously reported associations and DMRs identified in this study (Additional file 2: Table S2).The first two studies (Smith et al. (2021) [87] and Zhang et al. (2020) [86]) utilized methylation data from the Illumina HumanMethylation 450 k beadchip, from which we identified commonalities with our AD-associated DMRs (Additional file 3: Table S6).However, none of the 120 candidate ICRs overlapped the differentially methylated positions (DMPs) reported in these two studies.A limiting factor in using the Illumina HumanMethylation 450 k beadchip array is that coverage is for approximately 450,000 methylation sites, constituting only 3% of the total 28,084,558 CpGs in the human genome.Comparison of the 450 K manifest with the coordinates of The candidate ICR was previously defined by having 5 or more consecutive CpGs with methylation levels of 50% ± 15% (green dots) for tissues in all three germ layers (i.e., brain, kidney, and liver); methylation levels for sperm and oocytes are also shown (i.e., ≥ 90% methylation-yellow dots and ≤ 10% methylation-blue dots) [35] the AD-ICRs identified 102 CpG sites in common, in 45 ICRs.The lack of common regions between these different approaches is at least partially attributable to this low coverage.
The study by Breen et al. (2023) employed wholegenome methyl-sequencing in blood and looked at differentially methylated positions (DMPs) in AD patients compared to controls [88].The comparisons between the DMRs identified in this recent study and our AD associated DMRs listed in Table 1 revealed 102 overlapping DMRs in NHBs, including DMR chr7:76149513-76150704, which overlaps with candidate ICR_473, two DMRs in NHWs, and 22 DMRs in ALL group that are common in both studies (Additional file 3: Table S6).
Additionally, an epigenome-wide association study (EWAS) conducted by Piras et al. [89] revealed differential methylation patterns in the genes of AD brains as compared to non-demented control brains.The study identified 832 DMRs, out of which five DMRs were found associated with CASZ1, MAD1L1, MRPL23, RIMBP2, SLC9A3R2, ZCCHC14, and NFIC.Interestingly, these genes were each in proximity to an ICR identified in our study (ICR_20, ICR_434, ICR_716/ICR_719, ICR_805, ICR_922, ICR_976, and ICR_1071, respectively; see Table 1).There are other epigenome-wide studies that showed DMRs associated with various regions of the brain [87,90].However, more studies with a focus on ICRs and higher coverage are needed to understand the association between ICRs and AD development.
The AD-associated ICRs described here also overlap other known early-established methylation-dependent gene regulatory regions.A previous study determined regions of systemic interindividual variations (SIVs), characterized by consistent methylation across tissues within individuals, but with significant variation among individuals [91].Like ICRs, SIV methylation is established before tissue specification, hence the consistency within individuals, but unlike ICRs, methylation is not restricted to parent-of-origin status.The SIV control regions comprise approximately 0.1% of the human genome and regulate the expression of metastable epialleles.The most widely known example of a metastable epiallele is the agouti viable yellow (A vy ) locus in mice used to demonstrate that environmentally induced epigenetic modifications during early development can produce a range of phenotypes and alter disease susceptibility in adulthood [24].SIVs are defined as regions "conserved across diverse human ethnic groups, sensitive to periconceptional environmental exposures, and associated with genes implicated in a broad range of human disorders and phenotypes" [91,92].Interestingly, 15 of the AD-associated candidate ICRs overlap with previously described SIVs (Table 1) [91], including ICR_987 (chr17:5771207-5771575, NLRP1), common to both NHB-AD and NHW-AD case control comparisons.
Obesity is known as one of the modifiable risk factors for dementia [93].According to a study by Nianogo et al. [2] one third of the AD related dementia cases were linked to a combination of modifiable risk factors, including midlife obesity, physical inactivity, and educational attainment [2].In support of this, our analysis revealed the white adipose tissue browning pathway as a common enriched pathway associated with AD for both NHB and NHW populations.Changes in adipose tissue are part of the normal aging process [94], as adipose tissue is highly dynamic and has a role in homeostatic processes.White adipose tissue (WAT) may actively change into beige or brown adipose tissue (BAT) with environmental factors [95] through the white adipose tissue browning pathway.WAT is known to function as an endocrine organ secreting various types of adipokines including TNFα, which increases with aging.High WAT mass is related to metabolic disorders and linked to insulin resistance [94], and high BMI is associated with a reduction in brain volume [96,97].It is suggested that impairment in adipose tissue-derived adipokines may cause problems in brain homeostasis [94], which may eventually lead to neurodegenerative diseases.
Additionally, we identified the gap junction signaling pathway in both NHB and NHW populations as an enriched pathway, which is potentially dysregulated in AD.GNAS, an imprinted gene that is involved in multiple signaling pathways associated with G protein-coupled receptors, is also one of the molecules involved in the gap junction signaling pathway.Transcription of connexins, involved in the gap junction signaling pathway, are reported to be regulated by epigenetic modifications [98].GJA3, in close proximity to ICR_814 (chr13:20142811-20142911), is a member of the gap junction signaling pathway and was found to contain a differentially methylated CpG (chr13: 20736075) in AD hippocampus samples compared to controls [99].Future studies on the possible epigenetic regulation of these pathways may elucidate their mechanisms in the development of AD.
Interestingly, our network analysis identified netrin signaling in NHBs only.Netrins are axon guidance molecules associated with the regulation of axonal growth and play roles in neuroinflammation.Netrin 1 is reported to inhibit Aβ production [100].UNC5B, a molecule belonging to the netrin signaling pathway, was the closest gene to candidate ICR_664 (chr10:71266448-71266685), which was hypomethylated in AD brain samples compared to controls (Table 1, Additional file 2: Table S2).The UNC5B receptor is activated by the Netrin 1 molecule in the netrin signaling pathway.It is reported to have inhibitory roles in the inflammatory response in nervous system and may have a protective role in neurodegeneration [100,101].
A major strength of our study is that it is one the largest investigations, with deeply phenotyped participants with WGBS data, that uses DNA derived from brain samples not only from NHWs, but also NHBs, who are rarely included in the study cohorts, despite their higher risk of developing AD.Nevertheless, the data should be interpreted in the context of its limitations.Firstly, the relatively small sample size (n = 17) could diminish the statistical power required to identify some ICRs associated with AD.Despite the sample size, our use of agnostic WGBS revealed novel changes in ICRs that should be detectable in accessible tissues.Moreover, WGBS has revealed racial/ethnicity-dependent differences in ICR methylation that may lead to a better understanding of disparities in AD.Further analysis of this phenomenon may improve the treatment of AD by providing a mechanism for determining at-risk individuals.While the inclusive nature of our study aimed to encompass a diverse cohort, the small sample size also introduces an inevitable limitation, preventing us from drawing definitive conclusions regarding sex-specific ICR methylation patterns in the context of AD.However, further analyses with a larger sample size to address this issue in future studies are necessary.
Secondly, brain structural changes in AD start in the entorhinal cortex and medial temporal lobes and extend into the neocortex [4] and cerebellum [102] over time.We cannot exclude the possibility that the threefold difference in the number of differentially methylated ICRs associated with AD is due, in part, to epigenetic differences given the heterogeneity in NHB control tissues, which were a combination of temporal cortex and cerebellum.However, ICR methylation should be similar across tissue and cell types and should not be affected by tissue/cell heterogeneity.For that reason, we are confident with our definition of differentially methylation in candidate ICRs.
We also repeated our initial analyses using the more stringent DMR cutoff of 15% to reduce false positives and restricted our analysis to a subset of the 1488 ICRs previously reported, specifically the 332 ICRs with gametic methylation data.This more stringent analysis did not alter the findings that aberrant ICRs methylation is overrepresented among AD cases, more so in NHBs than NHWs.Although differential ICR methylation holds promise in surveillance to identify AD, replicating these findings in high-powered studies with DNA derived from various brain regions and accessible tissues such as blood or central nervous system fluids, in diverse populations, are necessary.
Finally, while we acknowledge the existence of epigenetic drift associated with aging, we lack information regarding the methylation status of ICRs throughout the aging process.ICRs represent specific regulatory regions, and the DNA methylation patterns established during the intrauterine stage remain conserved.Furthermore, a study conducted by Mancino et al. (2023) on the hippocampus of mice, highlights a notable age-related increase in DNA methylation-an established transcriptional indicator of aging [103].Interestingly, the authors observed that genomic imprinting, specifically parentof-origin-specific DNA methylation, remained largely unaffected by the aging process.This observed stability extended across various brain regions, including the cerebellum, nucleus accumbens, hypothalamus, and prefrontal cortex.Transcriptomic analysis further substantiated these findings, confirming the preservation of imprinted expression in the aged hippocampus [103].Nevertheless, our current dataset does not provide insights into agerelated changes.With the anticipation of acquiring more comprehensive data through larger sample sizes and long-term follow-up studies, we aim to unravel the methylation dynamics of ICRs in the future, drawing from an extended population.

Conclusion
Using unbiased WGBS, we provide the first evidence that DNA methylation in 120 ICRs varies markedly between AD cases and controls.The number of ICRs with altered methylation is three times higher in NHBs with AD than in NHWs with AD, which may contribute to the higher prevalence of AD in NHBs compared to NHWs [1].Our findings are also consistent with the developmental origins of health and disease (DOHaD) hypothesis that increased susceptibility to adult-onset chronic diseases such as AD frequently have their origins in early development [104], and support the findings that AD is characterized by changes in the brain that likely start decades before the clinical symptoms appear [40,41].Thus, alteration in ICR methylation may serve as an early detection tool of AD risk that is essential for slowing the progression of this disease.

Human brain specimens
To identify differentially methylated ICRs associated with AD, we obtained frozen autopsy brain specimens from the Joseph and Kathleen Bryan Brain Bank of the Duke University/University of North Carolina at Chapel Hill Alzheimer's Disease Research Center (Duke/UNC ADRC); informed consent documentation is in IRB ID# 00016278.The AD and control brain tissues were selected according to their neuropathologic diagnosis of AD.Nine brain samples from AD cadavers (five NHBs and four NHWs) and eight brain samples from control cadavers (four NHBs and four NHWs) were collected.For all individuals, the time elapsed between patient death and collection and snap-freezing of samples was less than 24 h (Additional file 1: Table S1).

Library preparation of specimens, WGBS, and identification of AD-associated DMRs
We performed WGBS to identify DMRs associated with AD in all brain samples (n = 17).Libraries were prepared from extracted DNA using EpiGnome ™ Methyl-Seq reagents (Illumina, Inc, San Diego, CA), index-tagged for multiplexing, and sequenced on an Illumina Next-Seq platform (Illumina, Inc, San Diego, CA).Reads were assigned back to individuals by index reads, and aligned in silico to a bisulfite-converted reference genome (i.e., hg38 version 87).Reads without unique alignments, due to either repetitive genomic sequence or loss of specificity from bisulfite conversion of cytosines, and duplicate reads, indicative of clonal amplification of original random DNA fragments, were eliminated.The quality of each paired-end sequence file was inspected using fastqc, and adapter trimming and quality control were performed using the Trim Galore wrapper script that calls the cutadapt [105] script internally.Paired-end reads with a trimmed adapter sequence were aligned to the human reference genome (i.e., hg38 version 87) downloaded from Ensembl using bsmap aligner [106].The following options were passed to the aligner: -p 8 -L 135 -w 100 -v 10 -q 10 -R -V 1.The aligned bam files were sorted and indexed using samtools [107], and duplicate sequence reads were removed using the Picard [108] application.
We analyzed three groups from 17 study participants whose clinic and demographic characteristics are summarized in Additional file 1: Table S1.These included: I) All AD samples vs. all controls, II) NHB-AD vs. Microsoft Excel v16.46 (21021202) and Bedtools v2.30.0 were used to compare ICR coordinates with published data.Using the Expression Database of Human Long noncoding RNAs (LncExpDB) [72], we compared lncRNAs associated with brain development genes, and the genes in close proximity to ICRs.Furthermore, using hsa.gff3data provided by miRbase [73], we identified microRNAs among the annotated genes that overlap with AD-associated ICRs.NHB-control samples, III) NHW-AD vs. NHW-control samples (Fig. 1a).Briefly, a methylation ratio for CpG dinucleotides was generated from replicate bam files for each AD samples and controls using mcall, and DMRs were called between AD samples and controls using mcomp.The following options were utilized for mcomp inputs: -doStrandSpecifiMeth = 1, -doDmrScan = 1, -doDmcScan = 1, -dmrMethods = 2 -minDmcsInDmr = 4, -minCredibleDif = 0.1 or 0.15, -maxDistConsDmcs = 300, -minDepthForComp = 7 -pFetDmc = 0.05.DMRs were called using Model-based analysis of bisulfite sequencing data (MOABS, version 1.3.8.7), which relies on Credible Methylation Difference (CDIF) as a single metric for both statistics and biological significance of differential methylation, i.e., significant DMRs were generated using MOABS, and confirmed with known DMRs as positive controls [42].The following criteria were used in calling a DMR: minimum differentially methylated C's ≥ 4, credible cis-acting methylation difference ≥ 10% or ≥ 15%, minimum read depth ≥ 7, and max distance between consecutive CpG's ≤ 300 (Fig. 4).After merg- ing the bam files, the total coverage from certain CpGs divided by total coverage for all CpG's (wig sum percentages) were plotted versus increasing read depth using mmint (https:// github.com/ lijia cd985/ Mmint) (Additional file 1: Fig. S1a).This plot helped us check if there were high duplication levels or sequence bias.There was no sequence bias observed for all AD and control samples.We have generated the coverage plot using the plot coverage function from deepTools [109] using the merged bam files per group (Additional file 1: Fig. S1b, c).

Identification of AD-related ICRs in NHBs and NHWs
The DMRs associated with AD were intersected with the recently defined ICRs [35] using bedtools.A description of the process for defining human ICRs was recently published [35].Briefly, puticr, a custom tool implemented in Python (Version ≥ 2.7), was used to identify 1488 ICRs in tissue from multiple germ layers [35].We used Bedtools [110] to identify ICRs that intersected with AD-associated DMRs, and then used two-sided Fisher's exact test to identity enrichment of ICRs in the set of AD DMRs (Fig. 4).
To explore the functional significance of these ADassociated ICRs, genes closest to either side of the 120 AD-associated ICRs were identified from the UCSC genome browser.To explore molecular and cellular functions, potentially enriched pathways, and associations with disease, we analyzed our gene set using online tools, including Ingenuity Pathway Analysis (IPA) [111].

Fig. 1
Fig. 1 AD-associated candidate ICRs in NHBs and NHWs.a DMRs that differed in DNA methylation (≥ 10%) between AD cases and controls in NHBs and NHWs were determined by WGBS.b Venn diagram of ICRs from ALL[40], NHB[81] and NHW[27] when DNA methylation differed by ≥ 10% between AD cases and controls.c Venn diagram of ICRs from ALL[10], NHB[32] and NHW[10] when DNA methylation differed by ≥ 15% between AD cases and controls.Created with BioRender.com

Fig. 3
Fig. 3 Race/ethnicity independent ICRs in AD. a ICR_481 (MEST/MESTIT1) and b ICR_987 (NLRP1) differed by ≥ 10% in DNA methylation between AD cases and controls in both NHBs and NHWs.Candidate ICR (horizontal red box) and a known ICR (horizontal yellow box) is delineated by vertical dashed red lines.The candidate ICR was previously defined by having 5 or more consecutive CpGs with methylation levels of 50% ± 15% (green dots) for tissues in all three germ layers (i.e., brain, kidney, and liver); methylation levels for sperm and oocytes are also shown (i.e., ≥ 90% methylation-yellow dots and ≤ 10% methylation-blue dots)[35]

Fig. 4
Fig. 4 Experimental workflow.The steps used to identify DMRs in AD cases vs controls, overlapping ICRs, and closest genes.Created with BioRender.com

Table 1
AD-associated DMRs overlapping with candidate ICRs